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The dynamical evolution of super star clusters (SSCs) moving in the background 
of a dark matter halo has been investigated as a possible event responsible for the 
formation of bulges in late-type spirals. The underlying physical processes include 
sinking of SSCs due to the dynamical friction and stripping of SSCs on their way 
to the center. Based on the assumption of a universal density profile for the dark 
matter halo, and an isothermal model for the SSCs, our simulations have yielded 
bulges that are similar in many aspects to the observational ones. In particular, 
the derived surface density profiles can be well fitted by an exponential structure 
with nuclear cusps, which is consistent with HST observations.The preliminary 
simulations with the Burkert density profile yield very instructive predictions on 
the bulge formation, which is surely worth further investigating. 


1 Introduction 

It is believed for a time that bulges in many late-type galaxies display exponential 
surface brightness profiles. Now the HST observations (Carollo & Stiavelli 1998; 
Carollo 1999) revealed a very interesting central structure in a samll sample of late- 
type galaxies, which illustrates typical surface brightness profiles of nucelar cusps 
on top of exponential bulges. Carollo and her colleagues convincingly claimed the 
stellar origin for these nuclear cusps. 

In fact, nuclear star clusters may be a common phenomenon (e.g. Davidge & 
Courteau 2002; and references therein) in late-type galaxies. Possibly related to 
the existence of nuclear star clusters is the presence of circumnuclear young massive 
star clusters, or circumnuclear SSCs. 

A question then arises whether there is a physical relation between these 
two common phenomena in late-type spirals, i.e., the common existence of nu- 
clear/circumnuclear star clusters and the exponential bulges with nuclear stellar 
cusps. If this is the case, what is the physical process responsible for this relation? 
A prevailing view for the formation of bulges in late-type galaxies is the secular 
evolution of disks due to the bar formation/destruction (e.g. Norman, Sellwood, & 
Hasan 1996). However, Bureau argued in his review (2002) against this scenario 
based on the consideration of the fast duty cycle of the bar formation/destruction, 
inferred from the omnipresence of bars. 

In this paper we report our study, motivated by Carollo’s claims and the ob¬ 
servations mentioned above, on the dynamical evolution of circumnuclear SSCs 
moving in the dark matter halo as a possible process for the formation of bulges 
with nuclear cusps in very late-type galaxies. 
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2 Models 


2.1 SSC model 

Based on the analogy with globular clusters (GCs), we assume that SSCs have a 
similar mass spectrum as the initial GC mass function but with larger mean value. 
We adopt the log-normal mass function as SSCs’ mass function (SSCMF) 

M 

Logio{——) ~ N{exp{mean) = 2 x 10 Mq, variance = 0.08) (1) 

Mq 

Considering the mean GC’s mass of about 1 — 2 x 10^Mq, we may infer the mean 
SSC’s mass of about 2 x 10®M©. 


2.2 Background 

Considering the bulge-formation process in very late-type galaxies (Sd type and 
later) we only take a dark matter (DM) halo as the initial background. For the 
dark halo, we assume the universal density profile (Navarro et al. 1997, hereafter 
NFW), which can be written as (e.g., Binney & Merrifield, 1998) 


Ph{r) 


Moh 

r{ah +rY 


( 2 ) 


where r is the distance from the halo center. 

We consider two halos with quite different masses: a large one with mass M 200 = 
lO^^ft.“^M 0 and a smaller halo with mass M 200 = 3 x IO^^Mq. The NFW halos 
are determined via the scaling law (see Table 1), with the Hubble constant Hq = 
75 kms“^Mpc“^, h = 0.75. 


Table 1. Parameters of the chosen NFW halos 


M20o{Mq) 

r 200 (kpc) 

c 

ah{kpc) 

Moh^Mq) 


216.8 

8.5 

25.5 

7.8 X 10^" 

3 X 10^® 

131.9 

10 

13.2 

1.6 X 10®" 


In order to account for the background variation caused by the fallen SSCs 
together with their stripped mass, we add a spherical component, represented by 
shells of equal density, which varies according to the changing SSC mass contribu¬ 
tion. 


2.3 Dynamical friction 

We define M and Vm (with Vm = \ Vm\) equal to, respectively, the mass and velocity 
of cluster experiencing the dynamical friction. Assuming the background matter 
has a Maxwellian velocity distribution with dispersion <Jbkgd, and is composed of 
particles with mass much smaller than M, the dynamical friction formula may be 
written as (e.g., Binney & Tremaine, 1987) 


dVM 

dt 



( 3 ) 
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where erf is the error function. The readers are advised on referring to Binney & 
Tremaine (1987) about A and X. 


2.4 Stripping 

We assume that the stellar mass outside a sphere of radius Rt, which corresponds 
to the instantaneous Hill stability region around the SSC center, will be stripped. 
As a first-order approximation, the mass of the stripped stars is considered to be, 
at some later epoch, uniformly distributed in the shell bounded by r -|- i?( (r) and 
r — dr — Rt{r — dr). By summing up all of the stellar mass stripped at various rs, 
the stripped stellar mass distribution can be derived. 

If a single massive object is embedded in the center of an SSC, stripping cannot 
proceeded further when only this object is left. Therefore, in our simulations, we 
consider two extreme cases: stripping is not allowed when the mass of the stripped 
SSC is less than IMq and IOOOMq, respectively. 


3 Results 

Figure 1 shows the mean surface density prohles of simulated galaxies with high- 
density central components formed at 1 Gyr (with a proportion of about 65% and 
90% for the larger and smaller halos, respectively). 

In order to give a quantitative comparison with the HST observations, the pre¬ 
dicted mean surface density profiles are fit with a model, following Carollo & Sti- 
avelli (1998), which is given by 


cr(i?) =croexp(-1.678^)-hcri(l-H^)'^exp(-^) (4) 

ite -^s 

The resulting prohles are also shown in Fig. 1, and the associated values of the 
parameters are given in Table 2. Except at radii less than about 2 pc, the httings 
shown in Fig. 1 are quite satisfactory. The half mass radii of the formed bulges are 
about 300pc, and the scale-sizes of the nuclear cusps are about lOpc , which are 
consistent with HST observations. 


Table 2. The htted values of parameters in & 


Af200 (Af©) 

(To {Mq/pc^) 

Re {pc) 

(Ti (MqIpc^) 

Rc {pc) 

7 

Rs {pc) 


127 

361 

317 

8 

2.2 

23 

3 X IQii 

146 

327 

212 

18 

1.7 

27 


In order to see what could happen to the galaxies where no high-density central 
component is formed at 1 Gyr, we extend the simulation to 1.5 Gyr. The bulge 
formation rates are then increased to about 90% and 99% for the larger and smaller 
halos, respectively. This result implies that the proposed mechanism may be a 
viable mode of bulge formation. 
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Figure 1. Mean surface densities derived by SSCs without IMBHs for two DM density profiles 
adopted. The model fittings are illustrated with solid blue and dashed red lines 


4 Discussion 

Now more and more evidence seems to argue (e.g. Marchesini et al. 2002 and the 
references therein) that the NFW profiles are not adequate for a fraction of dwarf 
galaxies which are dominated by the dark matter halos, and suggested that the 
Burkert profile is more suitable to these galaxies. We would then take the Burkert 
density profile given below for our simulations. 


PDM{r) 


Po r 


3 

O 


{r + ro){r^ + rl) 


( 5 ) 


The work of this kind has not been finished yet. One of our preliminary results 
is illustrated in Fig 2, obtained with a central density of 0.10 Mq pc~^ . It is very 
instructive to notice that the exponential bulges with nuclear cusps formed in about 
40% galaxies, and a large fraction of dwarf spirals have no bulges and no nuclear 
cusps formed. 
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Figure 2. Model fitting to the simulations with the Burkert density profiles 


Table 3. Comparisons of the predicted to the observed bulge formation rates 



NFW 

Burkert 

0.03 (po) 0.10 

Obs.(~Sd)* 

bulge 

size 

cusp 

size 

~90% 
300 pc 
~90% 
8—10 pc 

0% ~40% 

200pc 
0% ~40% 

12pc 

<10% - 50% 
N/A 

~30% - 50% 
<lpc — <20pc 


* based on Matthews et al. (2002, 1999) and Carollo (1999) 


In contrast to the case of NFW density profiles, the bulge formation rates yielded 
with the Burkert ones are obviously lower. To compare both of these predicted 
rates to the relevant data obtained for very late-type galaxies might give a clue 
to analysing the DM density profiles adequate to these galaxies. The comparisons 
summarized in Table 3 apparantly indicate that the Burkert density profile is more 
suitable to the galaxies observed by Matthews et al. (2002, 1999) and Carollo 
(1999). 


fhd: submitted to World Scientific on February 2, 2008 


5 














Nevertheless, we have to be cautious of the above claim. For one thing, the 
observed properties for galaxies with Hubble type of Sd or later are taken from 
a small sample of about 20 sources only. And more important is that the results 
obtained for the Burkert profiles are very preliminary. A detailed investigation on 
this topic is underway. 
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